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Asymptotic and chaotic solutions of a singularly perturbed 

Nagumo-type equation* 

Alberto Boscagginl Walter Dambrosio^ Duccio Papini^ 


Abstract 

We deal with the singularly perturbed Nagumo-type equation 

-I- m(1 — u){u — a{s)) = 0, 

where e > 0 is a real parameter and a : K —>■ K is a piecewise constant function 
satisfying 0 < a{s) < 1 for all s. We prove the existence of chaotic, homoclinic 
and heteroclinic solutions, when e is small enough. We use a dynamical systems ap¬ 
proach, based on the Stretching Along Paths method and on the Conley-Wazewski’s 
method. 
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1 Introduction 

In this paper, we deal with the singularly perturbed Nagumo-type equation 

e^u" + u{l — u){u — a{s)) = 0, (1-1) 

where e > 0 is a small parameter and o : M —>• M is a locally integrable function satisfying 

0 < a(s) < 1, for every s G M. 

Notice that equation (O) has the two constant solutions u = 0 and u = 1 which actually 
behave like saddle points); we will be interested in the existence of solutions u satisfying 
0 < u{s) < 1 for all s. 
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Our investigation is motivated by a classical paper by Angenent, Mallet-Paret and 
Peletier [T], dealing with the Neumann boundary value problem associated with 0 in 
the framework of steady-states solutions of the corresponding parabolic problem (which 
arises in population genetics). On the lines of previous works [H |T3], in [T] it is proved 


that, for e small enough, the Neumann boundary value problem associated with (1.1) has 
multiple solutions, whose shape and limit profile (for e —?■ 0"*“) can be precisely described 
in terms of the zeros of the function a — 1/2. This analysis suggests that the dynamics 


of (1.1) could be quite rich provided the function a crosses the value 1/2. 


The aim of the present paper is indeed to prove the existence of solutions to (1.1) 


defined on the whole real line and exhibiting complex behavior, when a switches inhnitely 
many times between two values a_, a+ with 

1 

0 < a- < - < a+ < 1. 


(see, however, Remark 4.5 for possible generalizations of the above assumption). More 


precisely, we will provide solutions of essentially two types. On one hand, we find the 
existence of globally defined solutions rotating (in the phase-plane) a certain number 
of times around (a±,0) in suitable intervals (j G Z) (solutions of this kind will be 


called chaotic)] the number of revolutions can be arbitrarily prescribed and Uj^zlf is 
unbounded from below and from above, i.e., solutions oscillate infinitely many times on 
the whole real line. In particular, if a is T-periodic for some T > 0, we find the existence 


of subharmonic solutions to (1.1) with complex behavior; indeed, one could show that 


the Poincare map associated with (1.1) on a period is topologically semiconjugated to 


the Bernoulli shift on a suitable number of symbols (as in [3], see also Remark 5.4). 

On the other hand, we are able to produce homoclinic and heteroclinic solutions, 
having the same nodal behavior as before on bounded intervals of arbitrarily large length 
and converging monotonically to one of the equilibria (0, 0) and (1,0) for t —)• Too. 

A similar nonlinearity is considered in the recent papers la [3 Eg EH]. However, 
those papers deal with situations which can be considered in a certain sense dual to 
ours, since they have a fixed central equilibrium (a center, in fact) while the two nearby 
saddle points move according to a suitable piecewise constant weight function. Moreover, 
we decided to focus on finding solutions that are obtained by exploiting quite different 
configurations than those investigated in the mentioned papers, and to omit the details 


in cases that are similar to those already treated (see Remark 4.6) 


For the proof of our results, we use the change of variable x{t) = u{et), so that the 


equivalent planar system associated with (1.1) is transformed into 


x' = y, y' = x{x — l)(x — a(et)); 


( 1 . 2 ) 
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hence, the smallness of e > 0 reflects into the fact that the function ae{t) = a{et) is 
constant on intervals of large amplitude. As a consequence, (1.2) can be regarded as 
a slowly varying perturbation of an autonomous system with two hyperbolic equilib¬ 
rium points; in this setting, the arising of rich dynamics seems to be a quite common 
phenomenum (see, among others, iHiiini as well as the bibliography in [29]). 


In order to detect such a complex behavior, we use a dynamical systems approach, 
based on a careful analysis of the trajectories of the (piecewise autonomous) system 
(1.2). More in detail, we first rely on a topological technique of path stretching (the 
so-called SAP method) developed in [^ EH ESI E3j to detect the presence of symbolic 
dynamics (as well as of periodic points, when a is periodic) for suitable Poincare maps 
associated with (1.2). As a by-product, this approach also provides us planar paths 
which can be used to connect the stable and unstable manifolds of the equilibria so as to 
obtain heteroclinic and homoclinic solutions with a complex nodal behavior (see, among 
others, PEI da IT5] and the references therein for related results in this direction). 
It has to be noticed that, since the equation is non-autonomous, the existence of the 
stable/unstable manifolds is not straightforward; we use indeed the classical Conley- 
Wazewski’s method mm (see also PIE]) to prove that these sets actually exist and 
can be suitably localized. 


The plan of the paper is the following. In Section]^ we briefly discuss the autonomous 
case. In Section]^ we prove some stretching properties, as well as the existence of stable 
and unstable manifolds to the equilibria. In Section [^ we state and prove our main 
results. Finally, some basic facts about SAP method are collected in a final Appendix. 


2 The autonomous case 

In this section we collect some basic results for the autonomous equation 

x” + x{l — x){x — a) = 0, (2.1) 

where a is a real constant such that 


0 < a < 1. 


Precisely, we are going to perform a phase-plane analysis for the equivalent planar system 


x' = y 

y' = x{x — l)(x — a). 


( 2 . 2 ) 


We will always confine our attention to the dynamics in the vertical strip [0,1] x M. 
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It is immediately seen that the points (0,0), (a, 0) and (1,0) are the only equilibrium 
points of (2.2). To proceed further, we define the function 


£a{x,y) =-y^ + Fa{x), V(x,y)G 


where 


Fa{x) = --x^ + 


1 T u 


-X ^x 


y X £ 


As well-known, system (2.2) is conservative and the function t i—)• £a{x{t),y{t)) is con¬ 
stant along solutions {x{t),y{t)) to (2.2). To describe the global dynamics of (2.2), we 
can thus study the geometry of the level sets for different values of c G M. 

It turns out that the value of the constant a plays a significant role. We start by 
analyzing the case 

‘ (2.3) 


“ 2 ’ 

which indeed gives rise to the simplest picture. Precisely, the level set can here 

be described as follows: 


- for c = ^^ 0 ?, it is the point (a, 0); 

- for ^^ 0 ? < c < 0, it is a closed cycle around (o, 0); 

- for c = 0, it is the union of the points (0,0), (1,0) and of the heteroclinic orbits 
joining them; 

- for c > 0, it is the union of two curves, one in the half-plane {(a:, y) G : y > 0}, 
one in the half-plane {(a:,y) G : y < 0}, connecting two points of the form 
(0,yi) and (l,yi) for some yi G M. 

The phase-portrait is shown in Figure 

We now turn our attention to the case 

0<a<^; (2.4) 

now, for the level set £~^{c) we have the following: 

- for c = ^^ 0 ?, it is the point (a, 0); 

- for ^^ 0 ? < c < 0, it is a closed cycle around (a, 0); 

- for c = 0, it is the union of the point (0,0) and its homoclinic orbit H[a)] for 
further convenience, we denote by {za,0) the point of intersection between H(a) 
and the positive x-semiaxis; 
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Figure 1: The phase-portrait of the autonomous system (2.2) for a = 1/2. The heteroclinic 
orbits connecting the equilibria (0,0) and (1,0) are painted with a darker color. For graphical 
reasons, a slightly different x and y scaling has been used. 


- for 0 < c < it is made by a curve lying between the homoclinic to (0,0) 

and the stable/unstable manifold H^[a) of (1,0), connecting a point of the form 
(0, yi) with a point of the form (0, —yi), for some yi > 0; 

- for c = it is the union of the point (1,0) and its stable/unstable manifolds; 

- for c > ^4^, it is the union of two curves, one in the half-plane {(x, y) G : y > 
0}, one in the half-plane {(x,y) G : y < 0}, connecting a point of the form 
(0, yi) with a point of the form (1, ^ 2 ), for some yi,y 2 G M. 

Finally, for 

^ < a < 1 (2.5) 

the phase-portrait can be obtained from the previous one with a symmetry with respect 
to the line x = ^. The homoclinic orbit and the stable/nnstable manifolds are defined 
in an analogous way, by swapping (0,0) and (1,0), and will be again denoted by H{a) 
and H^{a) ; moreover, ( 2 ;a, 0 ) will be the point of intersection between H{a) and the 
positive x-semiaxis. Both the phase-portraits are shown in Figure 
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Figure 2: On the left, the phase-portrait of the autonomous system (2.21 for 0<a<l/2;on 
the right, the phase-portrait of the autonomous system (2.2| for 1/2 < a < 1. The homoclinic 
orbits, as well as the stable/unstable manifolds, are painted with a darker color. For graphical 
reasons, a slightly different x and y scaling has been used. 


Remark 2.1. From the above discussion, it appears that the phase-portrait of system 
(2.2) is completely different in the case a = 1/2 and in the case a ^ 1/2. Indeed, the 
heteroclinc orbit connecting (0,0) and (1,0) for a = 1/2 disappear as soon as a / 1/2, 
splitting into orbits of different type. In terms of the potential Fa{x), we have indeed 
Fa{0) = F’a(l) if and only if a = 1/2; in this case, the potential is said to be balanced. In 
this context, (4.1) can be framed in the setting of equations with unbalanced potentials 
(compare with the introduction in 


3 Topological lemmas 

In this section we collect the preliminary technical lemmas which will be used in the 
proof of our main results. 


3.1 Stretching properties 

In this section, we fix two reals constants a_,a+ satisfying 


0 < a_ < - < a+ < 1. 


Our goal is to prove some results for the dynamics of the autonomous system (2.2) for 


a = a- and a = o+ on suitable sets which will be constructed below. Throughout this 
section, we always refer to the definitions given in the Appendix. Also, to simplify the 


notation, from now on we denote by S'(a-) (resp., 5'(a+)) the planar system (2.2) for 
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a = a- (resp., a = a+); moreover, let &{a,z) be the orbit of (2.2) passing through the 
point z G 

For T > 0 we finally define the maps and as the (restriction to the vertical 
strip [0,1] xM) of Poincare maps associated with systems S{a-) and S{a+), respectively, 
on the interval [0,T], i.e. 


^'±(xo,yo) = ix{T-,xo,yo),y{T;xo,yo)), V(xo,yo) G [0,1] x 


where {x{-]Xo,yo),y{-]Xo,yo)) is the unique solution to S{a±) satisfying the initial con¬ 
dition (x(0),y(0)) = {xo,yo). Since we will be interested in the dynamics on the strip 
[0,1] X M only, we can assume that such maps are globally defined just by suitably mod¬ 
ifying the nonlinearity fa{x) = x{l — x){x — a) for x ^ [0,1] (for instance, by setting 
fa{x) = 0 for X ^ [0,1])- 

Let us fix p- and p+ such that 

- < max {za_, a+} < < 1 and 0 < < min {a_, <-. (3.1) 

Let TZi and TZ^ be the two connected components of the intersection of the following 
two strips: 

- the strip S- between the stable manifold and the orbit 0(a_, (p_,0)) of 

5(a_), 

- the strip 5+ between the unstable manifold and the orbit 0(a_, (p+,0)) 

of 5'(a+). 

Notice that, since p+ < and Za_ < P-, the above defined orbits passing through 
{p±i 0) lie between the homoclinics and the stable/unstable manifolds of the correspond¬ 
ing systems S{a±). Therefore, the defined regions TZi and TZ^ are topological rectangles, 
and we name TZi (resp., TZs) the one contained in the upper (resp., lower) half-plane. 

In what follows we also need to provide an orientation for the above constructed 
rectangles; to this aim, we denote by TZf the components of the boundary of TZi lying 
on orbits of system S{a±) and by TZ^ the components of the boundary of TZi lying on 
orbits of system S{az^). 

Now, let 7^2 be the intersection of the following regions: 

- the strip S- defined above, 

- the annular region between the homoclinic H (a+) and the closed orbit 0(o+, (g+, 0)) 
of S{a+), for some q+ such that 


P- <q+ < 1 , 


(3.2) 
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- the upper half-plane. 


Straightforward computations show that the homoclinic H{a^) is contained in the region 
bounded by the stable and unstable manifolds Moreover, since a+ < p_ < 

(recall both (3.1) and (3.2)) the closed orbit 0(a+,(g+,O)) winds around the point 
(p_,0), as well. These facts together guarantee that 7^2 is a topological rectangle. 

In a similar way, we can define the rectangle TZ^ as the intersection of the following 
regions: 


- the strip defined above, 

- the annular region between the homoclinic H{a-) and the closed orbit 0(a_, (g_, 0)) 
of S{a-), for some q- such that 


0 < g- < P+, 


(3.3) 


- the lower half-plane. 


Notice that the rectangles 7^2 and TZ^ depend also on the choice of the numbers q± 


satisfying (3.2) and (3.3) which is not arbitrary and will be specified in Proposition 3.1 


We can again orientate the above constructed rectangles, denoting by 7^^ the com¬ 
ponents of the boundary of 7^2 lying on orbits of the systems 5(0+) and by 7^^ the 
remaining components; hnally, TZ~^ are the components of the boundary of 7^4 lying on 
orbits of the systems S{a-) and are the remaining components. 


El 


We illustrate the whole construction of the rectangles TZi, for z = 1,..., 4 in Figure 
We are now in position to prove our crucial result on the stretching: 


Proposition 3.1. The following stretching properties hold true. 


1 . There exists Tf > 0 such that, for every Ti > Tf, we have 

4/^^ : (7^l,7^^)^(7^2,7^2-), : (7^2,7^+) ^ (7^3,7^3-), 

: (7^3,7^3)^(7^4,7^J), : (7^4,7^+) ^ (7^l, 7^^), 

for a suitable choice of q'^ and q~ (close enough to 1 and to 0, respectively). 

2. For any iV G N, there exists T|(A^) > 0 such that, for every T 2 > T^iN), we have 

: (7^2,7^^)^'^(7^2,7^+), : (7^4,7^4 ) ^"^ (7^4,7^+). 
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Figure 3: The construction of the topological rectangles TZi, for i in particular, 

we comment in detail the construction of TZi and 7^2- Consider the stable manifold (to (1,0)) 
(in red), as well as the homoclinic (to (0,0)) H{a-) (in pink) and recall that such an 
orbit intersects the positive x-semiaxis at the point {za _, 0). The condition p_ > Za_ (see ( |3.1| )) 
then guarantees that the piece of orbit of 0(a_, (p_, 0)) (also in red) in the upper half-plane lies 
between H{a-) and H~^{a-). With an analogous construction for the system S{a+), and taking 
into account that p+ < 1/2 < p_ (see (3.1) again) we can thus determine the rectangle TZi (in 
orange) in the upper-half plane (as well as the rectangle TZ^ (in purple) in the lower half-plane. 
Now, consider the homoclinic (to (1,0)) H{a+) (in black), intersecting the positive x-semiaxis 
at the point {za^,0). It is easy to check that this orbit is contained in the region bounded 
the stable and unstable manifold moreover, the condition Za^ < a+ < p- (see ( |3.2[ )) 

implies that it intersects the orbit 0(a_, (p_,0)). Finally, focus on the orbit 0(a+, ((7+,0)) (also 
in black). Since a+ < p- < q+ (see both (3.1) and (3.2)), such an orbit intersects 0(a_, (p_,0)) 
as well (that is, the closed orbit 0(a+, (g+, 0)) winds around the point (p_, 0)). This determines 
the rectangle TZ 2 (in green). An analogous construction gives TZ 4 (in grey). Recall that, for the 
validity of the stretching properties in Proposition 3.1 q± cannot be arbitrary numbers satisfying 
(3.2) and (3.3), but they have to fulfill further conditions (see the proof). 
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Proof. 1. We give the details for the stretching property 

the other three statements of the first group being analogous. We first observe that 
the time needed to cover the piece of the orbit 0(a_,(p_,O)) contained in the right 
half-plane is given by 

_ . 

Jo y/Fa_{p-) - Fa_{xy 

similarly, the time needed to cover the piece of the orbit 0(a+, (p+, 0)) contained in the 
right half-plane is given by 


%/2 


dx 


'p+ 


Now, let 


Tf = V2 


max 





dx 


dx 


y/Fa_{p-) - Fa_ix)jp^ y/Fa+{p+) - K+(®) 


and fix Ti > Tf. Moreover, we define 

q+ = sup |x G [p_, 1] : 0) G 7^l| . 

Let us observe that the above set is non-empty since, in view of the choice of Ti, 
X I—)• 0) (x G [p_,l]) parameterizes a curve joining a point in the second 

quadrant with the point (1,0), whose intersection in the first quadrant is contained in 
S~. By a continuity argument, moreover, we can easily prove that p- < < 1. In this 


manner, (3.2) is satisfied and this completely determines the rectangle 7^2- This choice 


of helps in controlling the behavior of the solutions which move according to system 
S{a-) on a time interval of length Ti and either start from TZi or arrive on TZ^. In fact, 
such solutions evolve inside the strip S~ and, by the choice of qj^, they can cross the 
positive X axis only at a point lying between (p-,0) and ((jr+,0). Similarly, we define 

q- = inf |x G [0,p+] : (T+i)"^(x,0) G 7^3} . 
in order to complete the definition of 7^4. 


We now pass to the verification of the stretching property, according to Definition 5.1 
Let 7 : [0,1] —)■ 7^i be a path such that 7(0) and 7(1) lies on different components of 
TZf and, just to fix the ideas, suppose that 7(1) lies on the stable manifold 


10 












We observe that 'h^^( 7 (l)) remains on in the first quadrant while, by the choice 

of Ti, ^'^^( 7 ( 0 )) lies in the third quadrant; hence, by the positive invariance of the strip 
S- for the flow of 5(a_), we have 

^'!:l(7([O,l]))n7^2/0. 

On the other hand, by the choices of q+ and Ti, ( 7 ([ 0 ,1])) cannot intersect the 
boundary of 7^2 on the x-axis. As a consequence, it intersects the two components of 
7^2 ) say 7^2 ti '^2 d means “down” and u means “up”, with obvious meaning). 
Now, let 

t~ = sup |f G [ 0 , 1 ] : G 7^■^| 

and 

7+ = inf jf > 7“ : ^^^'( 7 ( 0 ) ^ ^ 2 ,«} • 

By construction, the sub-path 71 = o 7 l[t- t+] has image contained in TZ -2 and crosses 
it from 7^2 d ^2 ii> desired. 

2. We now turn our attention to the second group of statements, proving that 

: ( 7 ^ 2 , 7 ^ 2 -)^'^( 7 ^ 2 , 7 ^+) 

(the other one being analogous); we observe that here 7^2 appears with different orien¬ 
tations when considered as the domain or the target space of We introduce the 
polar coordinate system, centered at (a+, 0 ), 

X = a+-|-r cos 0 , y = —rsm6] 


then, for any 2 : G 7 ^ 2 , we denote by 0(-;z) the (unique) continuous angular function 
associated to the solution of S{a+) starting from 2 : at the time 7 = 0, and such that 
0{O]z) G [— 7 r, 0 ]. Moreover, we recall that the time needed by any solution to system 
5(a+) starting from a point on TZ^d perform one turn around (a+, 0 ) can be computed 
as 

^F,+ {q+)-F,+ {x) 

where G (0,g''') is the only number such that = F„+(g+). 

Given the natural number A7 > 1, we define 


T^{N) = {N + l)V2 [ 


dx 


4+ VFaAq+)-FaAx) 
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and let T 2 > T^. For any j = 1,..., we set 

nj = {z G 7^2 : 9{T2, z) € [{2j - l)7r, 2j7r]}; 
let us observe that the sets Hi,... ,'Hn are compact and disjoint. 

Now, let 7 : [0,1] —>■ 7^2 a path such that 7 ( 0 ) and 7 ( 1 ) lie on different components of 
7^2 and, to fix the ideas, assume that 7 ( 1 ) lies on 7^^.^ (that is, on the stable manifold 
to (1,0)). Hence, 0 (T 2 , 7 ( 1 )) < 0 while, in view of the choice of T 2 , 0 (r 2 , 7 (O)) > IN-n. 
As a consequence, we can find A'-disjoint subintervals /i,...,Ijv C [0,1] such that 
7 ( 5 ) G T-Lj for any s G Ij and o crosses 7^2 from one component of 7^^ to the 
other component of ^ 


Remark 3.2. We observe that the stretching relationships proved in point 2. of Propo¬ 
sition 3.1 naturally provide information about some nodal properties of the solutions 
of the systems S{a±). In particular, the (clockwise) winding number around the point 
(a+, 0) of solutions of S{a+) starting in Hj (for some j = 1, ..., A), in a time interval of 
length T 2 , lies in the interval {j — 1 / 2 , j -|- 1 / 2 ) and, more precisely, the first derivative 
of those solutions vanishes exactly 2j times. Similar considerations holds of course for 
solutions to S{a-) around the point (a_,0). With a slight abuse of terminology, in our 
main results we will say that, in this case, solutions make j turns around (a±, 0 ). 


Remark 3.3. We observe that the construction of the rectangles TZi (i = 1,...,4) 
given in this section is independent from the fact that the homoclinics H(a-) or F7(a+) 
intersect, that is, we are not assuming that Za+ < Za_ (see Figure]^. However, it is 
worth noticing that when the strict inequality 


^ ■2'a 


(3.4) 


holds true, one can easily construct two further rectangles satisfying stretching condi¬ 
tions. This is indeed a well-known geometrical configuration, related to the concept 
Linked Twist Map, which has been extensively discussed in m- We refer again to 
Figure]^ With elementary computations, it is possible to show that condition (3.4) is 
satisfied when a± are not too far from the value 1 / 2 . 


3.2 Stable and unstable manifolds 


In this section we prove the existence of stable and unstable manifolds to the equilibria 
of a nonautonomous planar system. Our result will be given in a slightly more general 
setting than needed; more precisely, we deal with the planar system 


x' = y 

y' = x(x - l)(x - g(t)), 


(3.5) 
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Figure 4: On the left, the construction of the rectangles when a_ = 0.4 and 0 + = 0.6, so 
that the homoclinics intersect (this is indeed the same picture as in Figure]^. On the right, 
the construction of the rectangles in a situation in which the homoclinics do not intersect; here 
a_ = 0.3 and a+ = 0.7. One can immediately realize that nothing changes. However, on the 
left again, it is shown how the linking between the homoclinics provides two further topological 
rectangles, painted in yellow. It is well-known that such rectangles verify stretching properties, 
within a Linked Twist Maps framework (see m)- 


with q : 


a locally integrable function. We use the notation z{t]p,to) for the 


(unique) solution to (3.5) satisfying the condition z{to) = p. 


Proposition 3.4. The following statements hold true. 


1. Assume that 


0 < a_ < q{t) < a+ < 1, for a.e. t G (— 00 , to]) 

for some to £ Then there exist two continua (i.e., connected compact sets) 
C lying between the unstable manifolds to (0,0) and (1,0) of the 
systems S{a-) and 5'(a+) respectively, and such that, for any p G r(_oo; holds 
that 

lim 2 ;(t;p, to) —)• (t, 0), for i = 0,1. 

t—)■—00 

Moreover, there exist two constants G (0, a_) and a\_ G (o+, 1) such that 
n ([0,a°] X M) = {(x,y°^(x)) : x G [0,a°]} 

and 

r]_^ n {[aX, 1] X M) = {(x, yi^ix)) : x G [a]^, 1]}, 
for suitable continuous functions y^ao ■ [0) ®-] ^-oo ■ [®+) 1] 


13 
















2. Assume that 


0 < a_ < q{t) < a+ < 1, for a.e. i G [to, +oo), 


for some to £ Then there exist two continua (i.e., connected compact sets) 
C lying between the stable manifolds to (0,0) and (1,0) of the 
systems S{a-) and S{a+) respectively, and such that, for any p G r^oo, holds 
that 

lim z{t-,p,to) ^ {i,0), fori = 0,1. 
t —>-+00 

Moreover, there exist two constants a® G (0,a_) and a\ G (a+, 1) (the same ones 
as in the previous statement) such that 

r+oo n ([0,a°] X M) = {(x,y[J.^(x)) : x G [0,o°]} 


and 

r^oo n {[a\, 1] X M) = {(x, ?/+oo(a;)) : a: G [a^, 1]}, 
for suitable eontinuous functions y%ao ■ [0; ®-] V+oo ■ [®+) 1] 


Proof. We observe that Statement 2 is a consequence of Statement 1 and the symmetry 
enjoyed by the vector field in (3.5) with respect to the axis y = 0. Hence we have that 
rVoo = {{x,y) ■ {x,-y) G r^oo} and y\^ = -yl^ for i = 0,1. Concerning the proof 
of Statement 1, we give only the details about the continuum since the existence 

and the properties of can be deduced in an analogous way. 

Let us fix e > 0 such that 


0<a_ — e<a+ + e<l 

and define C as the compact triangular region of the first quadrant of the phase- 
plane bounded by the unstable manifolds to (0, 0) for the systems S'(a_ — e) and S'(a++e) 
and the vertical line x = a_ — e. 

We first prove that any solution z{t) = {x{t),y(t)) which remains in for every t <tQ 
has to satisfy 

lim z(t) = (0,0). 

t^ — OO 

Indeed, for such a solution we have that, for t < to, 

x{t) G (0, a_ — e] , x{t) > 0, x"{t) > 0, 

and, hence, there exist 

lim xit) = L < a- and lim x'(t) = 0. 

t^ — OO t^ — OO 
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Finally, L must be zero; otherwise: 


liminf x”{t) = liminf{x(t)[l — x(t)][a(t) — x{t)]} > L{1 — a_)e > 0 

t^ — OO t^ — OO 

and, thus, we would have limi_,._oo x\t) / 0. 


Now we write system (3.5) as an autonomous system in 


X = y 

y' = x(x - l)(x - q{t)) 
t' = 1 


(3.6) 


and let 'k{-\sq,P) be the unique solution to (3.6) starting from P € at s = sq. We 
set 


VF = Se X (-00, to], U = (7- \ {( 0 , 0 )}) X (-00, to], V = (7+ \ {( 0 , 0 )}) x (-00, to]. 


where 7_ and 7+ are the portions of the boundary of lying on the unstable manifolds 
to (0, 0) for the systems S(a- — e) and S(a+ + e) respectively. 

Let us study the behavior of the vector field associated with system ( |3.6[ ) at any point 
P = (xi,yi,ti) G dW. First, if (xi,7/i) = (0,0), it is clear that 7r(s;ti,P) = (0,0, s) 
for all s. Next, if P G ?7 U 1/, the vector field points strictly inwards W] therefore, 
7r(s;ti,P) ^ W for all s in a left neighborhood of ti. Finally, if either 

(xi,yi) G \ (7_ U7+U{(0,0)}) 


or 

(xi,yi) G Se\ (7-U7 +U{( 0,0)}) and = 0 

then 7r(s;ti,P) G W for all s in a left neighborhood of ti, since in those points the 
vector field of (3.6) points strictly outwards W. 

Now, we consider the set D dW given by 

D = {P = {xi,yi,ti) G VF : 3s < s.t. 7r(s; p, P) ^ W} 

and the map 

D -^dW 


such that <h(P) is the first backward exit point from kF for the solution of (3.6) starting 
from the point P, namely <h(P) = '7r(s*;ti,P) where 

s* = sup{s < ti : 7r(s;ti,P) 0 D}. 
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It is proved in [5] that is continuous on D; moreover, by the previous discussion, 
$(D) = U L)V is not connected and U and V are its connected components. 

Let 7 : [0,1] —>■ be a continuous path such that 7(0) G 7_ \ {(0,0)} and 7(1) G 
7 + \ {(0,0)}. Then we have that <I>(7(0)) G U, <I>(7(1)) G V. Since [0,1] is connected, 
there must be r G (0,1) such that 7 (t) ^ D. By the topological lemma [25l Corollary 
6] there exists a continuum T*}C \ ZD such that 

(0,0) G and n ({o_ - e} x M) ^ 0. 

Letting e —)■ O"'', it is actually possible to show that C So (that is, it lies between 
the unstable manifolds to (0,0) and (1,0) of the systems S{a-) and 5(a+) respectively) 
and reaches the vertical line x = a_ (see [T3l Theorem 6, §47, II, p.l71]). 

The fact that, in a suitable vertical strip, r[{;^ can be written as the graph of a continuous 
function can be proved as in [26( Lemma 2.4], taking into account the fact that, for 


f{t,x) = x{l — x){x — q{t)) and a^_ = 
it holds dxf{t,x) > 0 for a.e. t < to and x G [0,a*}]. 


1 T ct— — \ ci_ — U— ~\~ 1 


□ 


Remark 3.5. We observe that the mere existence of the stable and unstable manifolds 
follows from the Stable Manifold Theorem fsee|lljL since (0, 0) and (1,0) are hyperbolic 


equilibrium points of system (3.5). Here, we have used Wazewski’s method in order to 


provide, for such sets, a precise localization, which is indeed needed in our arguments. 
Notice that in this way we cannot directly claim that stable/unstable manifolds are 
curves (indeed, Wazewski’s method applies in more general cases in which this is not 
true); however, we can recover such an information in an elementary way, using the sign 
of the nonlinearity as in 


4 Main results 

Let us consider the equation 

e^u" + /(s, u) = 0 

where 

f m(1 - m)(u - a(s)) if'uG[0,1], 

if [0.1], 

with a : M —>■ M a locally integrable function satisfying 

0 < a(s) < 1, for every s G M. 


(4.1) 

(4.2) 
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Remark 4.1. We observe that, in view of the boundedness of /, any solution to (4.1) 
is globally defined. Moreover, if a solution u satisfies u{s*) ^ (0,1) for some s* G M, 
then u{s) ^ (0,1) either for all s < s* or for all s > s*. 


From now on, we assume 

a(s) = 

where a_, a+ are real constants with 


(X— if S2k — S ^ ^2/c+l 

«+ if S 2 k +1 < s < S 2 k +2 


0<a_<-<a+<l 


and {sk)k&z C M is a sequence such that: 


0 < (5 < inf {sk+i - Sk) 


(4.3) 


(4.4) 


(4.5) 


and, thus, Sk < Sfc+i for every /c G Z and limfc_,.±oo Sk = ±00. For every j G Z let us set 
k ~ [^6(7 —1)) ^6^]) Ij — ['S6j— 5 ) 'S 6 j_ 4 ], Ij = [sgj—2)'S6j—l]- 
In this setting, we will prove the following results. 


The first one. Theorem 4.2, deals with chaotic solutions. 

Theorem 4.2. For any integer M > 1, there exists e* = e*{M) > 0, such that, for any 
e G (0, e*) and for any double infinite sequence n = nj)}ji=z, 'with G {1,..., M} 
for all j G Z, there exists a globally defined solution u of ( |4.1[ ), with 0 < rt(s) < 1 for 
all s G M, such that its trajectory {u,u') makes turns around (a±,0) in the interval 
ifj for oil j G Z. Moreover, if the sequence {sk — Sfc-i) is 6-periodic and the sequence 
n is £-periodic for some I" G N, then this solution u is £(se — sq)- periodic. 

Notice that the sentence “n^ turns around (a±,0)” in the above statement has to 


be meant according to Remark 3.2; the same terminology will be used in Theorems 4.3 
and 4.4 below, dealing with heteroclinic and homoclinic solutions, respectively. 

Theorem 4.3. For any integer M > 1, there exists e* = e*{M) > 0, such that, for any 
e G (0,e*), for any integer K > 0 and for any n = {{n'^,n~)}j=i^,,,^K (this K-uple is 
considered to be empty if K = Oj with G {1, ..., M} for all j = 1,..., k, there exists 
a globally defined solution u of (4.1), with 0 < u{s) < 1 for all s G 

1. (n(— 00), tt'(— 00)) = (0, 0), 

2. (u(+cx)),u'(+cx))) = (1,0), 


such that: 
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3. the trajectory {u,u') makes turns around (a±,0) in the interval , for all 


4- u is monotone in (—00, sq] o,nd in [s6i^,+cx)). 


Theorem 4.4. For any integer M > 1, there exists e* = e*{M) > 0, such that, for any 
e G (0,e*), for any integer K > 0 and for any n = {(n^, (this K-uple is 


considered to be empty if K = Oj with G {1,..., M} for all j = 1,.. 


a globally defined solution u of (4.1), with 0 < u{s) < 1 for all s G 


k, there exists 
such that: 


1. (u(—cx)), m'(— cx))) = (0, 0), 

2. (u(+cx)),m'(+cx))) = (0,0), 


3. the trajectory {u,u') makes turns around (a±,0) in the interval , for all 


4- u is monotone in (—cx),so] and in [s 6 A'+i,+oo) and u' vanishes exactly once in 
(s6K, S6K+i)- 


4.1 Proof of the main results 

In order to prove the above theorems, we perform the change of variable 


x(t) = u(et), t G M, 


converting the equation (4.1) into 


x" + f(et,x) = 0, 


with 


f(et,x) = < 


x(l — x)(x — o_) if X G [0,1] and t G t 2 k+i) , /c G Z 
x(l - x){x - a+) if X G [0,1] and t G [t 2 fc+i,^ 2 ^+ 2 ) ,k €Z 
0 if X 0 [0,1] 


where we have set 




Sk 


for every A: G Z. 


Notice that, in view of (4.5), we now have 


0 < - < inf (4+1 - 4). 

6 k£Z 
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Let us define, for any k € Z 


T^k = 'Pk = Til and Vk = ) 


and 


_ ^*6fe+6-i6fe+5 Q ^t6k+5-t6k+4 ^ ^j,*6fc+4-t6fe+3 ^ ^j^t6fe+3-*6^+2 ^ 


O ^*6*+2“*6*;+l Q ^hk+l-tek 


In other words, ^k{p) = {x{tQk+6), x'(tek+e)) where x is the solution of x” + f{et, x) = 0 
with {x{tek),x^{hk)) =P- 

Proof of Theorem \4-^ Given M > 1, let us define e*{M) > 0 as 


e*{M) = 


max{rf,r|(M)}’ 


(4.6) 


where Tf ,Tf{M) are given by Proposition 3.1 Now we take e G (0,e*); notice that we 
have 

tk+i — tk > maxlT^, r2 (M)}, for every k G Z. 


In view of Propositions (3.1) and ( |5.2[ ), we have 

4 : Vk Vk- 


The conclusion thus follows from Theorem 5.3 in the Appendix. We finally notice that 
x{tQk) G (0,1) for all k G Z, so that, in view of Remark 4.1, x{t) G (0,1) for all t G M. □ 


Proof of Theorem \4-3[ We begin by applying Proposition |3 . 4| on the intervals (—oo, t_i] 
and [t6fc+i,+oo) in order to find the continuous functions ■ [Oj®-] and 

y+oo : [«+>l] such that (x,y°^(x)) G and {x,yl^{x)) G P^^oo for all 

suitable x. We denote by x* the abscissa of the intersection point in the first quadrant 
between the homoclinic to (0,0) for the system 5'(a_) and the orbit 0(a+,(p+,O)) of 
system 5(a^). Let 

xi = min{x*, a®}. 


The time needed by a solution of system S{a+) to run along 0(o+, (p+, 0)) from x = xi 
to X = 1 is given by 

1 dx 

\/2 Jxi •\/Va+ {p+) ~ Va+ ( 3 ^) 

In a similar way, we take x** as the abscissa of the intersection point in the first quadrant 
between the homoclinic to (1,0) for the system S{a+) and the orbit 0(a_,(p“,O)) of 
system S{a~). Let 

r ** It 

0:2 = maxjx 
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The time needed by a solution of system 5(o_) to go on 0(a_, {p ,0)) from x = 0 to 
X = X2 is given by 

, _ 1 _dx_ 

\/2 Jo y/Fa_{p-) - Fa_{x) 

For any integer M > 1, we can thus define 


\M) = 


maxlTj^, r2* (M),r, r'}’ 


(4.7) 


where again and T^iM) are given by Proposition 3.1, With these positions, we have 
that 






(r_oo 0 ([0>a:i] X 


is a path which crosses connecting the two components of TZ^ . 

Indeed, this is a consequence of the fact that (0,0) is a fixed point for and that 

the point (xi,yi) G F^^^ is mapped by in the half-plane {x > 1} which follows 

from the choice of r and the monotonicity of the map 


c I-)- 



dx 

\/c-Fa+{x) 


Analogously (by the choice of r') 


‘6'=)-i(r^^n([o,xi] xM)) 


is a path which crosses connecting the two components of . 

Now, corresponding to the choice of n = {(n^,with G {1,...,M}, 
there exists a sub-path 7 of *“^(F°g^) which is stretched across (7^i,7^^) by the 
map (t)k o 4>k-i o 

By the topological lemma m Lemma 3] the intersection 

(4 o 4_y o... o 4)(7) n 


is not empty and any point in this intersection gives rise to the required solution. □ 

Proof of Theorem 4-4 can argue in a similar way as in the proof of Theorem 4.3 
using here the curves F® and given in Lemma 


3.4 


□ 


Remark 4.5. We notice that the statements of our results can be modified (in some 
cases) to cover also more general equations e^u” + u{l — u){u — d{t)) = 0 with d close 


to a stepwise function. For instance, in the setting of Theorems |4.3| and 4.4 and given 
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integers M, K and e G (0, e*{M)), the existence result still holds for all functions a with 
L^-norm on [0, sqk] smaller than a constant 6 = 5{M, K, e) (compare with [161 Remark 
4.1] for more details on the stability of the stretching technique, and recall that Lemma 


3.4 about the existence of stable/unstable manifolds is indeed proved for a more general, 
non-stepwise, function). We also observe that the value e*{M) is an explicit constant 


(see (4.6) and (4.7)). 


Remark 4.6. We collect here some hints for possible variants and generalizations of 
our results. 


The role of the equilibria (0,0) and (1,0) may be switched in Theorems 4.3 and 


4.4 


More precisely, we can provide also solutions homoclinic to (1,0) as well as 
heteroclinic solutions converging to (1,0) for t —)■ —oo and to (0,0) for t —)• +oo. 
Moreover, the arguments used in the proof of Theorem |4 . 3| can be easily modified in 
order to find multiple solutions for Sturm-Liouville like boundary values problems 
(e.g., Dirichlet and Neumann ones). 


According to the discussion in Remark 3.3, it is possible to produce a further 


family of chaotic solutions when condition (3.4) is satisfied. 


In our main results we have considered functions of the form (4.3) with a_ and 


a+ satisfying condition (4.4), suggested by the investigation in [T]. When either 
a+ < 1/2 or a_ > 1/2, the superposition of the phase-portraits of the systems 
5(0+) and S{a-) gives rise to a different configuration (see Figure]^ analogous 
to the one considered in the paper [29], where a Schrodinger equation is studied. 


As a consequence, by combining the arguments therein together with Lemma 3.4 


it is possible to obtain the existence of chaotic dynamics and homoclinic orbits. 


5 Appendix: SAP method 

In this appendix we collect the definitions and results on the Stretching Along Paths 
method which are needed in our paper. We refer to |2l |22| for a comprehensive presen¬ 
tation of the theory and further references. 

By a path 7 in we mean a continuous mapping 7 : [0,1] —)■ M^, while by a sub¬ 
path cr of 7 we just mean the restriction of 7 to a compact subinterval of [0, Ij. By an 
oriented rectangle we mean a pair 

TZ = [TZ, TZ ), 

being 7^ C homeomorphic to [0,1]^ (namely, a topological rectangle) and 

TZ = TZ^ U 7^2 
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Figure 5: The superposition of the phase-portraits of the systems S'(a_) and S{a+), with 
a_ = 0.3 and a+ = 0.4. Orbits of the system S'(a_) are painted in pink, orbits of the system 
5(0+) are painted in blue, topological rectangles verifying stretching properties are painted in 
orange. Compare with Figure 2 in [2^ . 


the disjoint union of two compact arcs (by definition, a compact arc is a homeomorphic 
image of [0,1]) TZ^,71^ C dTZ. 

With these preliminaries, we can give the following definition. 

Definition 5.1. Let A = {A,A~), B = (B, B~) he oriented rectangles and let 'f' : T>^ C 
—>■ be a continuous map. 

- We say that ('H,'!') stretches A to B along the paths and write 

:A^B 

if 71 C .An Vq, is a compact subset and for every path 7 : [0,1] —)> .4. such 
that 7(0) G Af and 7(1) G Af (or 7(0) G Af and 7(1) G Af ), there exists a 
subinterval [t',t”] C [0,1] such that for every t G [t',7'] 

Jit) G n, ^ilit)) G B, 

and, moreover, 'I'(7(t')) and ^ij{t")) belong to different components of B~. 
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- We say that 'I' stretches A to B along the paths with crossing number M > 1 and 
write 

^ : A^^B 

if there exist M pairwise disjoint compact sets Hi,... ,'Hm C ^ n such that 
{Hi, ; A^B fori = l,...,M. 

As an easy consequence of the definition, we have that the stretching property has 
a good behavior with respect to compositions of maps. 

Proposition 5.2. Assume that 'h stretches A to B with crossing number M and 
stretches B to C with crossing number N. Then, the composition <ho 'I' stretches A to C 
with crossing number M x N. 

We finally state the result which is employed in our paper. 

Theorem 5.3. Assume that there are double sequences of oriented rectangles {Vk ) 

^ ^ V /fcez 

and of maps {4>k)k£Z that fk stretches Vk to Vk+i with crossing number Mk > 1, 

for all k G Z. Let Hkj C Vk, with j = 1,... ,Mk, be the compact sets according to the 

definition of multiple stretching. Then, the following conclusion hold: 

• for every sequence {sk)k&z, with G {1,..., Mk}, there exists a sequence {wk)kez 
with Wk G ^~Lk,s^, and (jkiwk) = Wk+i for all k G Z; 

• if there exists h,k G Z with h < k such that Vh = Vk, then there is a finite sequence 

{wi)h<i<k with Wi G Hi,Si, (fiiwi) = Wi+i for i = h,... ,k - 1 and Wh = Wk, that 
is, Wh is a fixed point of 4>k-i ° ° 4>h- 

Proof. By assumption, we have that 

{Hk,SkAk) ■ Vk^Vk+i, for every k G Z. 


The conclusion follows then from |21l Theorem 2.2]. 


□ 


Remark 5.4. In the setting of Theorem 5.3, when Vk = V, (fk 


k G Z, it is possible to prove that there is a compact invariant 
map 01A is topologically semiconjugate on the Bernoulli shift 
Lemma 2.3]). 


= 0 and Mk = M for all 
set A C P such that the 
on M symbols (see [221 


23 



References 


[1] S.B. Angenent, J. Mallet-Paret and L.A. Peletier, Stable transition layers in a semilinear boundary 
value problem, J. Differential Equations 67 (1987), 212-242. 

[2] L. Burra, Chaotic Dynamics in Nonlinear Theory, Springer, 2014. 

[3] A. Capietto, W. Dambrosio and D. Papini, Superlinear indefinite equations on the real line and 
chaotic dynamics, J. Differential Eqnations 181 (2002), 419-438. 

[4] Ph. Clement and L.A. Peletier, On a nonlinear eigenvalue problem occurring in population genetics, 
Proc. Roy. Soc. Edinburgh Sect. A 100 (1985), 85-101. 

[5] C. Conley, An application of Wazewski’s method to a non-linear boundary value problem which 
arises in population genetics, J. Math. Biol. 2 (1975), 241-249. 

[6] E. Ellero and F. Zanolin, Homoclinic and heteroclinic solutions for a class of second-order non- 
autonomous ordinary differential equations: multiplicity results for stepwise potentials. Bound. 
Value Probl. 2013, 2013:167, 23 pp. 

[7] E. Ellero and F. Zanolin, Connected branches of initial points for asymptotic solutions, with appli¬ 
cations, in preparation. 

[8] P.L. Felmer, S. Martinez and K. Tanaka, High-frequency chaotic solutions for a slowly varying 
dynamical system, Ergodic Theory Dynam. Systems 26 (2006), 379-407. 

[9] A. Gavioli, Monotone heteroclinic solutions to non-autonomous equations via phase plane analysis, 
NoDEA Nonlinear Differential Equations Appl. 18 (2011), 79-100. 

[10] T. Gedeon, H. Kokubu and K. Mischaikow, Chaotic solutions in slowly varying perturbations of 
Hamiltonian systems with applications to shallow water sloshing, J. Dynam. Differential Equations 
14 (2002), 63-84. 

[11] J.K. Hale, Ordinary differential equations. Second edition. Robert E. Krieger Publishing Co., Inc., 
Huntington, New York., 1980. 

[12] P.J. Holmes and C.A. Stuart, Homoclinic orbits for eventually autonomous planar flows, Z. Angew. 
Math. Phys. 43 (1992), 598-625. 

[13] K. Kuratowski, Topology, Vol. 2, Academic Press, New York, 1968. 

[14] H.L. Kurland, Monotone and oscillatory equilibrium solutions of a problem arising in population 
genetics, in: Nonlinear partial differential equations (Durham, N.H.,1982), Contemp. Math. 17, 
323-342, Amer. Math. Soc., Providence, R.I., 1983. 

[15] A. Margheri, C. Rebelo and F. Zanolin, Connected branches of initial points for asymptotic BVPs, 
with application to heteroclinic and homoclinic solutions, Adv. Nonlinear Stud. 9 (2009), 95-135. 

[16] A. Margheri, C. Rebelo and F. Zanolin, Chaos in periodically perturbed planar Hamiltonian systems 
using linked twisted maps, J. Differential Equations 249, (2010), 3233-3257. 

[17] J.S. Muldowney and D. Willett, An elementary proof of the existence of solutions to second order 
nonlinear boundary value problems, SIAM J. Math. Anal. 5 (1974), 701-707. 


24 



[18] K. Nakashima and K. Tanaka, Clustering layers and boundary layers in spatially inhomogeneous 
phase transition problems, Ann. Inst. H. Poincare Anal. Non Lineaire 20 (2003), 107-143. 

[19] D. Papini and F. Zanolin, A topological approach to superlinear indefinite boundary value problems, 
Topol. Methods Nonlinear Anal. 15 (2000), 203-233. 

[20] D. Papini and F. Zanolin, Fixed points, periodic points, and coin-tossing sequences for mappings 
defined on two-dimensional cells, Fixed Point Theory Appl. 2004, 113-134. 

[21] D. Papini and F. Zanolin, On the periodic boundary value problem and chaotic-like dynamics for 
nonlinear Hill’s equations, Adv. Nonlinear Stud. 4 (2004), 71-91. 

[22] A. Pascoletti, M. Pireddu and F. Zanolin, Multiple periodic solutions and complex dynamics for 
second order ODEs via linked twist maps. The 8th Colloquium on the Qualitative Theory of 
Differential Equations, No. 14, 32 pp., Proc. Colloq. Qual. Theory Differ. Equ., 8, Electron. J. 
Qual. Theory Differ. Equ., Szeged, 2008. 

[23] M. Pireddu and F. Zanolin, Fixed points for dissipative-repulsive systems and topological dynamics 
of mappings defined on N-dimensional cells, Adv. Nonlinear Stud. 5 (2005), 411-440. 

[24] M. Pireddu and F. Zanolin, Cutting surfaces and applications to periodic points and chaotic-like 
dynamics, Topol. Methods Nonlinear Anal. 30 (2007), 279-319. 

[25] C. Rebelo and F. Zanolin, On the existence and multiplicity of branches of nodal solutions for a 
class of parameter-dependent Sturm-Liouville problems via the shooting map, Differential Integral 
Equations 13 (2000), 1473-1502. 

[26] A.J. Urefia, Invariant manifolds around equilibria of Newtonian equations: some pathological ex¬ 
amples, J. Differential Equations 249 (2010), 366-391. 

[27] T. Wazewski, Sur un principe topologique de I’examen de Failure asymptotique des integrales des 
equations differentielles ordinaires, Ann. Soc. Polon. Math. 20, (1947) 279-313 (1948). 

[28] C. Zanini and F. Zanolin, An Example of Chaos for a Cubic Nonlinear Schrodinger with Periodic 
Inhomogeneous Nonlinearity, Adv. Nonlinear Stud. 12 (2012), 481-499. 

[29] C. Zanini and F. Zanolin, Complex dynamics in one-dimentional nonlinear Schrodinger equations 
with stepwise potential, preprint (2014). 


25 



